A critique of q-entropy for thermal statistics 
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During the past dozen years there have been numerous articles on a relation between entropy 
and probability which is non-additive and has a parameter q that depends on the nature of the 
thermodynamic system under consideration. For q = 1 this relation corresponds to the Boltzmann- 
Gibbs entropy, but for other values of q it is claimed that it leads to a formalism which is consistent 
with the laws of thermodynamics. However, it is shown here that the joint entropy for systems having 
different values of q is not defined in this formalism, and consequently fundamental thermodynamic 
concepts such as temperature and heat exchange cannot be considered for such systems. Moreover, 
for q 7^ 1 the probability distribution for weakly interacting systems does not factor into the product 
of the probability distribution for the separate systems, leading to spurious correlations and other 
unphysical consequences, e.g. non-extensive energy, that have been ignored in various applications 
given in the literature. 
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Introduction 

In 1988 a relation between entropy and probability for 
thermal statistics was proposed by C. Tsallis |l]] which is 
non-additive and depends on a parameter q that is pre- 
sumably determined by the nature of the thermodynamic 
system under consideration. For the special case q = 1 
this relation reduces to the Boltzmann-Gibbs entropy, 
but for other values of q it is claimed that this q-entropy 
leads to an alternative formalism that is "entirely consis- 
tent" with the laws of thermodynamics pi ||, [| ||, ||, |) . 
It will be shown, however, that this claim is not valid, be- 
cause the total entropy of thermodynamic systems that 
are in thermal contact is not even defined in this for- 
malism for the case where systems have q-entropies with 
different values of q. As a consequence, the fundamen- 
tal concepts of temperature and heat exchange between 
such systems can not be introduced in this formulation 
of thermodynamics unless q is a universal constant. But 
even in this case, maximizing the q-entropy of the com- 
bined systems leads to unphysical results unless q=l. In 
particular, for q =/= I the resulting joint probability for 
states of weakly coupled systems is not the product of 
the individual probabilities, and as a consequence the to- 
tal energy of the combined systems is not the sum of 
the mean energies of each system, although the additiv- 
ity property is assumed to be satisfied by the individual 
micro-states. This non-additive property of the mean en- 
ergy in the q-entropy formalism is manifestly incorrect, 
and not surprisingly leads to unphysical results, but the 
consequences have been ignored in various applications of 
q-entropy in the literature, as will be illustrated in some 
of the examples discussed here. 

The definition of q-entropy for a thermodynamic sys- 



tem with micro-states labelled by an index i is given by 

ill! 



bq ~ k (g-l) ' 



(1) 



where k is a constant, q is an undetermined parameter, 
and the quantities pi are positive numbers which satisfy 
the condition YliPi = 111 the limit that q — + 1 one 
recovers the Boltzmann-Gibbs form of the entropy 
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where k = ks is the Boltzmann constant. In this special 
case, pi is the probability for the occurrence of the i-th 
micro-state, and this identification has been extended to 
the case g ^ 1. This extension, however, is not valid 
as can be seen from the definition of mean values for 
physical quantities associated with the q-entropy. For 
example, the internal energy U q is given by the form 
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(3) 



where ej corresponds to the i-th energy eigenvalue of the 
system, and Pj =p\l Y^j P% The quantities Pi are called 
"escort" probabilities in the q-entropy literature, but ac- 
cording to the conventional definition of mean value in 
statistics, these P/s are the actual probabilities for the 
states of the system. Hence for q ^ 1 the p\s introduced 
in the definition of q-entropy, Eq. [I], are devoid of any 
physical meaning, and are just functions of the probabil- 
ities Pj according to the relation 
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For values of q ^ 1, the q-entropy expression intro- 
duced in Eq. |l| is shown to be non-additive by the fol- 
lowing arguments pl-p} . Suppose that two thermody- 
namic systems A and B are weakly coupled or are the 
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subsystems of a larger system, and assume that the joint 
probabilities for the states of the combined system are the 
product of the probabilities for the states of the individ- 
ual systems. As we shall see later on, this fundamental 
factorization property is not satisfied by the q-entropy 
formalism, but surprisingly this fact has been ignored in 
the literature. According to Eh. we then have 
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(5) 



Substituting this form into the expression for q-entropy, 
Eq. [l], one obtains the relation 

S q (AB) = S q (A) + S q (B) + (1 - q)S q (A)S q (B)/k. (6) 

However, this non-additive relation for the q-entropy 
leads immediately to a difficulty in the interpretation of 
S q as an expression for the thermodynamic entropy 0. 
Since weak coupling means that the energy eigenvalues 
of the combined systems are essentially additive, we have 



e AB 



(7) 



and, according to Eqs|| and |5[ the total mean energy 
U q (AB) of the combined system is also additive: 



U q (AB) = U q (A) + U q (B) 



(8) 



Assume now that the combined system is isolated while 
there is an infinitesimal exchange of energy between sys- 
tems A and B. Then the variations 5S q (AB) — and 
5U q (AB) = 0, which implies that 



6S q (A) 



SS q (B) 



1 + (1 - q)S q (A)/k l + (l-q)S q (B)/k' 

and 

SU q (A) = -5U q (B). 
Combining these two equations, one finds that 



(9) 



(10) 



[1 + (1 - q)S q (A)/k]T(A) = [1 + (1 - q)S q (B)/k}T(Bm) 

where T(A) and T(B) are the absolute temperatures de- 
fined by the standard thermodynamic relation 
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dU 
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(12) 



For two systems in thermal contact these two tempera- 
tures should be equal, but according to Eq.[y] this ap- 
plies only if q = 1. This problem is not unexpected, 
because temperature should be an an intensive quantity, 
but this is not possible in a formalism where the energy 
is additive, Eq.|lO|, while the entropy does not satisfy this 
property. To avoid this problem, it has been proposed 
||, ||, [h]] to re-define absolute temperature as the quan- 
tity 



in which case the condition for thermal equilibrium, 
T q (A) = T q (B) , is satisfied by Eq. 0. But this defi- 
nition of temperature, which must be universal, cannot 
be extended to systems A and B which are described by 
q-entropies with different values qA and qs ■ In this case 
the q-entropy of the combined system, which is charac- 
terized by the quantities pff, Eq. ||, is not defined in 
term of these two q parameters. Following the q-entropy 
formalism, one would have to introduce a new parame- 
ter q' for the the q-entropy and energy of the combined 
system, but then these thermodynamic variables cannot 
be expanded in terms of the corresponding variables for 
the component systems A and B ftl] , p2] |. Instead, ac- 
cording to Eq|| the q-entropy of the combined system 
would be given in terms of pseudo q-entropies for systems 
A and B with the same new parameter q' , and a simi- 
lar problem would occur with the expansion of the total 
energy, Eq.[| i.e. 

S q ' {AB) = S q/ (A) + S q ,(B) - (1 - q')S q ,(A)S q ,(B)/jm) 
and 



U q i{AB) = U q '(A) + U q i(B). 



(15) 



Consequently, the concepts of thermal equilibrium, tem- 
perature, and heat exchange cannot be formulated for 
such systems. For example, no meaning can be at- 
tached to the statement that a system described by the 
Boltzmann-Gibbs entropy is in thermal equilibrium with 
a system described by q-entropy with q =/= 1. In other 
words, a Boltzmann-Gibbs thermometer would not be 
able to measure the temperature of a q-entropic system, 
and the laws of thermodynamics would therefore fail to 
have general validity jl3| . 

It follows that the parameter q must be a universal 
constant, just like the Boltzmann constant k, which is 
applicable to all systems in thermodynamic equilibrium. 
If q is universal, a thermodynamic formulation for an 
infinitesimal reversible transfer of heat dQ can be given 
between systems A and B, with 



dQ = T(A)dS q (A) = -T(B)dS q {B), 



(16) 



corresponding to an exchange dU q (A) = —dU q (B) in 
the internal energy of these systems. But one is faced 
with the problem that in this case the temperature 
T(A) is not equal to T(B), which violates a fundamen- 
tal principle of thermodynamics for systems in thermal 
equilibrium. Moreover, the corresponding differentials 
T q (A)dS q (A) and T q (B)dS q (B) associated with the pro- 
posed re-definition of absolute temperature, Eq.|l3|, do 
not have any physical significance. In principle, this 
problem can be solved p"2] ] by introducing a different form 
for the q-entropy in the thermodynamic relation for 
the temperature T q , so that 



T q = (1 + (1 - q)S q /k)T 



(13) 
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From the definition of T q given by Eq. [L3| it follows that 
k 



oR _ 
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+ (1 - q)S q ), 



(18) 



This alternative expression for the q-entropy was intro- 
duced by A. Reyni [0| in the form 



(1-9) 



(19) 



which is additive |44j, as can be verified by substituting 
for the p^s the form in Eq. pi 

However, there are additional problems with either of 
these two definitions of q-entropy, even when qA = Qb = 
q. Maximizing such entropy functions subject to the con- 
straint of constant energy, Eq. || yields the q-probability 
distribution @, |, |, f§ 



(20) 



where Sei = €i — U q , U q is the mean energy, Eq. [| an d /3 is 
a parameter related to the inverse temperature. It is clear 
that the corresponding distribution pf^ for the combined 

system AB does not factor into the product pfpf even 
when the energies of the micro-states are additive, Eq. 
|?], unless q = 1. For example, to first order in q — 1 Ej| 
one finds that 

p(ef + ef) = p(e?)p(ef)[l + (q- l)(3 2 (5ef5ef)] (21) 

In the limit q = 1 Eq. po| reduces to the Boltzmann-Gibbs 
exponential form for the probability MM, 



Pi oc exp(-eJkT), 



(22) 



where (3 = 1/kT. As is well known, this canonical dis- 
tribution follows uniquely from the factorization require- 
ment that 



p(e?)p(ef)=p(e? + ef). 



(23) 



The Boltzmann-Gibbs form for the entropy S, Eq. |^, 
then follows from the assumption that 



(24) 



where the function f(p) is determined uniquely. Ap- 
plying the thermodynamic definition of temperature, 
dS/dU = 1/T and Boltzmann-Gibbs definition of the 
probability distribution, Eq. E2l one obtains the relation 



/ 1 1 



Hence 



dpi 



f{pi) = -kpilnipt 



(25) 



(26) 



provided /(l) = /(0) = 0, which corresponds to the re- 
quirement that the entropy vanishes at T = 0. 



Unphysical properties resulting from applications of 
q-entropy to thermodynamic systems 



To illustrate the consequences of disregarding such 
basic considerations, I would like to call attention to 
some unphysical results, left unmentioned in the liter- 
ature, that follow from recent applications of q-entropy 
to some well-known thermodynamic systems. For ex- 
ample, many papers have been published on the ^appli- 
cation of q-cntropy to black-body radiation jl5| [Ig[ IT], 
[0| |l9[ |(| |lj |||. After laborious analysis and numer- 
ical computations, the authors in references |l^, ^(| ||l[] 
find that for 9^1 there are deviations from the well- 
known Stefan-Boltzmann law which states that the radi- 
ation energy depends on the fourth power of the temper- 
ature jl6| . But since Boltzmann derived this result from 
purely thermodynamical reasoning, without any statisti- 
cal assumptions about the form of the entropy, it seems 
at first sight strange that such a deviation can occur in a 
formalism which is supposed to satisfy the laws of ther- 
modynamics. The explanation is that as a consequence 
of the non-extensivity property of the q-entropy formal- 
ism, the black-body energy as well as the entropy do not 
dependent linearly on the volume of the cavity, as was 
originally assumed by Boltzmann. Hence, the q-energy 
density and the q-entropy density depend on the volume 
V of the cavity, although there is no comment concerning 
this unphysical property in any of the q-entropy calcu- 



19], [|], H |§. As 



lations in references |15J, |16|, |17|, |18[ 
follows from simple dimensional arguments, the volume 
dependence must be given by a dimensionlcss parame- 
ter V/a(T) 3 , where a(T) is a characteristic length which 
can depend only on the temperature T. From statistical 
mechanics we learn that a(T) = (hc/kT) in fact cor- 
responds to the mean thermal wavelength of the black- 
body photons. Moreover, expansions in a power series 
of this parameter which have been applied to fit the cos- 
mic background radiation (l| [jj], y, [l|, §(| || || 
are nonsensical, because in this case a{T) is of order 1/10 
cm while the cavity volume V has cosmological dimen- 
sions! On purely thermodynamic grounds it can also be 
shown that if the temperature dependence of the black- 
body energy density were to have the form u oc T 4+s , 
then the Maxwell-Boltzmann relation p = (1/3)m for 
the thermal radiation pressure p would lead to a power 
law volume dependence u oc V s and correspondingly 
U oc T(T 3 V) 1+S/3 , in accordance with our previous di- 
mensional argument. Similarly, one finds that s oc V s / 3 
for the entropy density. But unless S = 0, such a volume 
dependence is incompatible with Kirchhoff's law which 
states that the ratio of emissivity to absorption of radia- 
tion in the walls of the cavity must be a universal function 
of the temperature, and the frequency of the radiation. 
These properties are required in order that such a cavity 
reach thermal equilibrium. Historically, this law was the 
original basis for the universality properties of black-body 
radiation that culminated in Planck's famous derivation. 

Actually, the application of the q-entropy formalism 
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to any system which consists of non-interacting compo- 
nents leads to unphysical properties. For example, for an 
ideal gas the energy is found to depend nonlinearly on the 
number of particles, the gas pressure p is not equal to 2/3 
of the energy density, and there are correlations between 
the energies of any pair of particles |23|, |2~i| , 25 1, contrary 
to very well-known result of kinetic gas theory and statis- 
tical mechanics. Similarly, for a system of non- interacting 
magnetic moments anomalous dependences are found for 
the magnetization and susceptibility on the number of 
spins [^6| , although claims have been also made that the 
magnetization reproduces experimental results in certain 
manganites |27], The reason for the failure of this for- 
malism to give physically sensible results is not hard to 
see - it is due to the fact that the probability distribu- 
tion, Eq. |2^, for any two components A\ and A2 is not 
the product of the probabilities for the separate com- 
ponents as would be expected if these components are 
non-interacting. Hence, contrary to basic physical prin- 
ciples, this formalism gives rise to spurious correlations 
among these components of such systems. 

Some proponents of q-entropy have argued that this 
formalism should be considered only for systems for 
which the Boltzmann-Gibbs thermodynamic formalism 
supposedly "fails". Frequently mentioned as candidates 
are systems with components that interact primarily 
through long range forces such as gravitational forces for 
which the total energy and entropy are non-extensive. In 
Nature, such system correspond to astrophysical objects 
such as stars and galaxies, but these objects are gen- 
erally not found at maximum entropy and (correspond- 
ingly) uniform temperature. For example, stars are ei- 
ther evolving slowly in time, like main sequence stars, 
which have a large temperature gradient from the inte- 
rior to the surface, and emit thermal radiation (black- 
body energy and entropy), or else have reached a degen- 
erate state such as dwarf stars or neutron stars provided 
certain mass limits arc satisfied |2j|. Otherwise, stars 
eventually explode into supernovas, sometimes leaving 
remnants which collapse into such degenerate states or 
into a black hole (which are states having maximum en- 
tropy). The equilibrium properties of stars are obtained 
by hydrostatic equations supplemented by local thermo- 
dynamic equations for matter and radiation based on the 
Boltzmann-Gibbs entropy |H| [^?J Rather than being a 
"failure", the Boltzmann-Gibbs statistics has been ap- 
plied to stellar structure with enormous success. A q- 
cntropy formalism, however, predicts the existence of fi- 
nite isothermal polytropes as the end products of stellar 
evolution ( states of maximum entropy) which fails com- 
pletely to account for the observed property of stars in 
our universe. Such q-polytrope solutions have also been 
discussed as models for galaxies 32, [33l B4J. 



of the hydrogen atom 



In this case 



Z = ^2exp(-e n /T), 



(27) 



Another example which has been cited as a so-called 
failure of Boltzmann-Gibbs entropy because it involves 
long-range electromagnetic forces, is the divergence of 
the partition function Z calculated for the bound states 



where e„ = —R/n 2 are the bound-state energy levels, 
and R is the Rydberg constant. But the divergence of 
Z , which occurs in this case because the terms of the 
series approach unity as n becomes large, is related to 
the growth of the the mean radius of the hydrogen atom 
which increases as n 2 . Obviously, in a gas of hydrogen 
atoms in thermal equilibrium, this radius can not become 
larger than the mean distance between atoms. Therefore 
this distance provides an effective cutoff for the applica- 
bility of the hydrogen bound-state energy eigenvalues in 
the partition sum, because for larger values of n these 
atoms can no longer be treated even approximately as a 
gas of non-interacting particles. Instead, for these states 
the gas must be viewed as a neutral plasma of electrons 
and protons interacting via long range electromagnetic 
forces. Hence the q-entropy formalism, which supposedly 
gives a finite partition function p5| , actually fails to ac- 
count for the correctphysics of this problem. These and 
other failures []36f , ]4q| . in the application of q-entropy 
to well-known physical systems mirror the inconsistencies 
which are inherent in a formulation of thermodynamics 
based on q-entropy, Eq. [l| 



A. Conclusion 

We have shown that a prerequisite to have a q-entropy 
formalism which is consistent with the laws of thermo- 
dynamics is that the parameter q must be a universal 
constant, as is the case also with the constant k that cor- 
responds to Boltzmann's constant for q = 1. Moreover, 
for weakly coupled systems the entropy as well as the 
energy must be additive, which is a condition explicitly 
violated by the Tsallis entropy, Eq. |]. Although this 
condition appears to be satisfied by the Reyni q-entropy, 
Eq. |l9], this is actually not the case, because the probabil- 
ity distribution obtained by maximizing this entropy for 
q 7^ 1, Eq. |2^, does not satisfy the required factorization 
condition, Eq. ^ which is required for weakly coupled 
systems. Thus, we have shown that the only value of 
q consistent with the laws of thermodynamics is q = 1, 
which corresponds to the familiar Boltzmann-Gibbs form 
for the entropy, Eq. ^. Indeed, it has been explicitly 
demonstrated here that the application of a q-entropy 
formalism to black-body radiation, and to other systems 
with weakly interacting components, leads to unphysical 
results when q ^ 1. 

It has been suggested in the literature that q-entropy 
calculations are useful because they provide an additional 
parameter q for comparing theory with observations, but 
this rationale fails to take into account the fact that such 
calculations would be inconsistent with the fundamental 
principles of thermodynamics and statistical mechanics. 
Hence, if a small departure of q from unity is found in a fit 
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to data, as is claimed for example in various analyses of 
the cosmic black-body radiation @, ||, ||, g|, |§ ||, 
EH], |2^], such a fit can not provide any physical insight 
whatsoever into the source or meaning of the deviations. 

Finally, we remark that the q-probability distribution, 
Eq. 0' which is obtained by maximizing either of the 
two q-entropy relations, Eqs. [I] and |l^, constrained to 
a fixed pseudo-energy function, has also been applied to 
non- equilibrium problems. For example, data on turbu- 
lence in a pure-electron plasma column [ |37| and the ve- 
locity distribution in turbulent flow |38| |3SH have been 
fitted by such a probability function 1 4C , [^J . But abso- 
lutely no physical justification has been given for apply- 
ing to non-equilibrium systems a fundamental condition 
- maximum entropy - which is associated in statistical 
mechanics with systems in thermal equilibrium. While 



such a probability distribution can also be obtained from 
other ad-hoc assumptions po| , |4lf| , its connection to q- 
entropy and non-extensivity is completely unfounded. 
Although the q-probability distribution appears to be a 
good phenomenological parameterization for some tur- 
bulence data, its deduction from physical principles has 
not been established . 
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cently, however, the SNO data (Q.R. Ahmad et. al., 
Physical Review Letters 89 (2002)) show that there is 
now excellent agreement between theory and neutrino ob- 
servations. Among other things, this result is a spectac- 
ular new confirmation of the applicability of Boltzmann- 
Gibbs entropy in the interior of a star 

[49] In references [^, the author expressed the q- 

probability distribution, Eq. as a Laplace integral, 
which is the basis for his suggestion that this distri- 
bution is due to temperature fluctuations in ordinary 
Boltzmann-Gibbs statistics. But no connection based on 
physics has been made between these assumed tempera- 
ture fluctuations and q-entropy 



